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^ ' Abstract. In this paper, radiating fluids scaling laws are studied. We focus on optically thin 

and optically thick regimes which are relevant for both astrophysics and laboratory experiments. 
I ■ By using homothetic Lie groups, we obtain the scaling laws, the similarity properties and the 

P ' number of free parameters which allow to rescale experiments in the two astrophyscial situations. 

ti ' 

1. Introduction 

' High-Energy-Density Physics is a new way for astrophysicists to explore phenomena usually 

l/^ ! occurring in the Universe. The use of powerful facilities, enables us to bring the matter up 

' to extreme states of density and temperature in laboratory [l]. The astrophysical relevance of 

. these experiments can be checked from scaling laws provided the physical system under study 

*0 I satisfies similarity properties. Thus, scaling laws and similarity properties must be examined 

■ with rigorous formalism. Several studies have been published about similarity and scaling laws. 
, For instance, in |2], [3], purely hydrodynamics and MHD scaling laws are respectively considered 

and in [2], the Birkhoff polytropic symmetries |4j are recovered. Moreover, in optically thin 
radiative hydrodynamic scaling laws have been considered and in [6j, the author not only studied 
similarity in case of optically thin plasma too but, also, discussed non-LTE situations through 

■ a microscopic approach. All these works have been carried out in an astrophysical context 
and were mainly based upon dimensional arguments. Scaling laws were also obtained for the 
Inertial Confinement Fusion (ICF) [7] , [8] in order to determine the minimum energy required for 
ignition. These are very interesting too because they can be used as non trivial tests for numerical 
simulations. In this paper, we study the radiating fluid similarity problem in two different 
regimes that can be (or will be) achieved in laboratory with current or future facilities. In each 
case, we derive the corresponding scaling laws and in order to get rigorous and exact relations, 
our approach is based on the Lie groups [9]. In the first part, we describe this method and remind 
its fundamental concepts. The second part deals with the optically thin radiating fluids, which 
are a major topic in astrophysics. Comparisons with other results obtained earlier are carried 
out. Finally we consider the equilibrium diffusion approximation including radiative pressure 
and energy. For each approximation, connections with astrophysical objects are provided and 
we emphasize the number of free parameters left to rescale an experiment. 



2. Lie groups, similarity and scaling laws 

The invariant transformation group theory elaborated by Sophus Lie is a very powerful tool of 
theoritical physics to study the symmetry properties of partial differential equations (PDE) and 
to perform their analytical integration. Among all Lie groups, one of them, namely, the one- 
parameter homothetic group (HG) is frequently used, first because of its simplicity and, then, 
because it provides more general self-similar solutions than those derived from dimensional 
analysis. This property arises because the HG is a sub-group of scaling transformations. 
Now, remembering the philosophy of Laboratory Astrophysics (i.e. to recreate systems having 
astronomical size on short scales), it seems natural to use the HG in order to study similarity 
properties, scaling laws and even self-similarity. Here, we will focus on the first two points 
only. Group invariance of PDE together with their solutions implies that the initial conditions 
(IC) be preserved from the laboratory system to the astrophysical one. This intuitive but 
constraining condition is discussed in details in [2], [3j. From this property, experimental data 
and IC in laboratory provide, first, information about astrophysical environments and, second, 
a transposition to astronomical objects. Moreover, the invariance of equations by HG implies 
that the Rankine-Hugoniot relations are also invariant and, therefore, we make sure that small 
scale shocks correspond to the homothetic structures of astrophysical shocks. 



3. Similarity and scaling laws of optically thin radiating fluids 

When the cooling (or heating) characteristic time of a plasma gets close to its dynamical time, 
this should be considered in the modeling. Concerning optically thin plasmas, i.e. Xp >> L (Ap 
is the mean free path of photons and L is the characteristic plasmas of the system), a simple 
modeling of radiating losses (or heating) can be done simply by introducing a loss (or gain) of 
entropy. Thus, the plasma is described by the following equations: 

^ + V.[pv]=0, P^ = -^PtH, = -(7-1)^/^'^)' dM = p.dV, (1) 

where d/dt is the Lagrangian derivative and p, v, Pth, 7 and M are respectively the density, 
velocity, thermal pressure, polytropic index and the mass of the fiuid. The function C{p, T) 
writes C{p,T) = Qi(p, T) + Q2{p,T) where Qi and Q2 are energy sources (or losses). In 
addition, we assumed a polytropic evolution; i.e. Pth = (7 — l)pe where e is the specific 
internal energy. Finally, an equation of state should be added to close ([T]): Pt^ = eo[Z]p>^T'^ 
where eo[Z] is a function of the ionization Z. It should be noticed that to satisfy the first 
thermodynamical principle we should have 7(1 — ly) = {p — v). Experimentally, heating can 
represent the laser energy deposition. From an astrophysical viewpoint, this modeling describes 
interstellar jets, bow shocks, radiating shocks (point C in Drake diagram [10], Fig 7.17) in Polars 
and supernova remnants. The relation between the typical quantities in astrophysical objects 
and laboratory experiments (that we note with ~) are given by: r = a'^if, t = a^^i, v = a^'^^v, 
M = a^^M, p = a^'p, Pth = a^'Pth, Qi = a^'Qi, Q2 = a^'Q2, eo = a^'io, T = a^^of, 7 = a'5"7 
where a is the group parameter and 6i are the homothetic exponants. Rescaling eq, and Qi 
can absorb a modification of Z from one system to the other (for example in bremsstrahlung 
cooling Q (X Z'^ ). Up to now, the sources have not been specified but in the applications 
we will consider power law forms {Qi = Qo,j/o'"'i'"'r''). This type of source is quite suitable 
for cooling since several processes in the continuum write in this simple form. We can also 
write Qi oc KpaT'^ where a is the Stefan-Boltzmann constant and Kp is the Planck mean 
opacity that can be modeled by a power law at high temperature. The invariance of equations 
under the HG provides the group invariants [9] namely: Ii = vt/r = St (Strouhal number), 
I2 = 7, Is = Ptht/pvr = Eu X St = St/ljM"^] {Eu\ Euler number, M: Mach number), 
Ii = Qit/Pth oc t/tQj^, Is = Q2t/Pth oc t/tQ^, where tg- is the characteristic time of the sources 
Qi and Ig = M/[pr^~^'^] (mass conservation). As expected, the invariants of this group are 



Table 1. Scaling for optically thin plasmas for power law models of sources (second column). Plane 
(d—0) radiative shock problem for magnetic cataclysmic variables: the third column corresponds to 
Bremsstrahlung Cooling (BC) [A oc p^T^''^] which can be Chevalier- Imamura unstable Ulf and the fourth 
column is obtained for BC plus cyclotronic cooling ( CC) [A oc p and a = PthP ■ 
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identical to the dimensionless numbers derived in similarity studies [6j. However, our approach 
is more general since we have local dimensionless quantities in contrast to global dimensionless 
numbers obtained thanks to the dimensional analysis. Thus, in our extension, the physical fields 
are conserved. Table 1 shows scaling laws for polars. Generally, we have four free parameters 
{5i, (55, ^6 and 5q) and, if we preserve ionization, only two (resp. one) exponent(s) remain(s) for 
a single source (resp. two sources). Moreover, if we set 6^ = 0, 6i = 1, 6q = 2 and Q2 = 0, we 
get the similarity considerations of [12] . Thus, with the same formalism, we can study similarity 
properties, scaling laws, and include the specific case presented in [12]. 

4. Similarity and scaling laws of optically thick radiating fluids 

Many systems, as well in laboratory as in astrophysics, are optically thick to radiation. For 
instance, the many classes of stars are more or less affected by radiation. Radiative pressure 
implies that there is an upper limit to the mass of a star (Eddington limit). Generally, including 
radiative fiux in laboratory experiments is enough and that is why, researches about scaling laws 
in this regime have been carried out in ICF. Here, we add the energy and pressure of radiation ( 
see [13] ) in the diffusion approximation at ETL. In experiments, LTE is usually satisfied [6j and 
it will be achieved on LMJ and NIF. In Astrophysics, radiation pressure and energy play a key 
role in stars, supernovae, in evaporation phenomena, in clumps |14]... The plasma evolution is 
then governed by the equations ([10], pp 270-271): 

dv -t , D 1 d f , TP \ pe + Pth + Erad + Prad dp - - n fO\ 

P-^ = -^[Pth + Pradh -^[pe + Erad) = -V.F^ad " Q, (2) 

where F^ad, Erad, Prad and Q are respectively the radiative fiux, radiative energy density, 
radiative pressure and the energy source term. In the application, we will consider that 
Erad = clrT^, Prad = Erad/ 3, Erad = -i^rad^T where Krad is the radiative conductivity given 
by Krad = Ko/o'"T"[we Still have Pth = Sqp^^T^]. In addition to the optically thin case we add 

the radiative relations: Erad = a^^^Erad] ^rad = a^^'-^^rad] = a^^^ko; Erad = a^^^ Erad] Prad = 
a^^'^Prad', Q = a'^^^Q. As before Ii, I2, I3, I4 (or I^) and Iq are five invariants. The additional 



Table 2. Scaling laws of optically thick plasma (Column 1) and ideal gas (Column 2). 
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ones are I'j = Pradt/ipvr) = Eurad ^ St {Eurad- Radiative Euler number), /g = Erad/Pth ocl/R 
(R: Mihalas numbers), /g = tFrad/ {Pth''') = I/-B0 {Bo: Boltzmann number). We recover the 
standard dimensionless numbers [13] which describe these radiating fluids. The scaling laws are 
presented in table 2. We see that for an ideal gas with an ionization conservation state, we have 
a single (^5) parameter to rescale experiments, but if the ionization is not preserved, we have 
at least three free parameters (^5, Jg, (5i4). Finally, we find that a (entropy) is conserved for 
7 = 4/3, which corresponds to a dominant photon regime. Notice that if we set Prad = Erad = 0, 
we find an extended scaling laws version of [7]. 

5. Conclusion 

We presented the general scaling laws for two radiative regimes of major interest in laboratory 
and astrophysics situations. The number of free parameters depends on the structure of the 
model: the more phenomena we add, the more difficult it is to rescale an experiment. However, 
requiring a partial similarity {^almost' equivalent regime) allows to add free parameters and 
study 'almost' astrophysical situations. 
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